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Background: Good understanding of the cross sections for (anti)neutrino scattering off nuclear targets in the 
few-GeV energy region is a prerequisite for the correct interpretation of results of ongoing and planned oscillation 
experiments. 

Purpose: Clarify a possible source of disagreement between recent measurements of the cross sections on carbon. 

Method: Nuclear effects in (anti)neutrino scattering off carbon nucleus are described using the spectral function 
approach. The effect of two- and multinucleon final states is accounted for by applying an effective value of the 
axial mass, fixed to 1.23 GeV. Neutral-current elastic (NCE) and charged-current quasielastic (CCQE) processes 
are treated on equal footing. 

Results: The differential and total cross sections for the energy ranging from a few hundreds of MeV to 100 GeV 
are obtained and compared to the available data from the BNL E734, MiniBooNE, and NOMAD experiments. 

Conclusions: Nuclear effects in NCE and CCQE scattering seem to be very similar. Within the spectral function 
approach, the axial mass from the shape analysis of the MiniBooNE data is in good agreement with the results 
reported by the BNL E734 and NOMAD Collaborations. However, the combined analysis of the NCE and CCQE 
data does not seem to support the contribution of multinucleon final states being large enough to explain the 
normalization of the MiniBooNE-reported cross sections. 

PACS numbers: 13.15.+g, 25.30.Pt 



I. INTRODUCTION 

We are currently witnessing a very rapid progress 
in the determination of neutrino oscillation parameters. 
The mixing angles and the mass-squared differences gov- 
erning oscillations of solar [lHH and atmospheric 0, HI 
neutrinos are known with systematically improving pre- 
cision; for a recent review see Ref. @. 

The past-year results from the T2K 0, MINOS @, 
and Double Chooz [Sj experiments have brought indi- 
cations that, in the three-neutrino framework, the last 
unknown mixing angle is nonvanishing. These hints have 
been confirmed by the findings from Daya Bay [loj . fol- 
lowed by those from RENO [11], with statistical signifi- 
cance of 5.2 and 4.9 standard deviations, respectively. 

Three nonvanishing mixing angles are a prerequisite 
for the existence of the phase violating the charge-parity 
(CP) symmetry. Therefore, the results reported by the 
Daya Bay and RENO Collaborations are paving the way 
for studies of CP violation in the lepton sector, which, in 
turn, may bring us closer to understanding the matter- 
antimatter asymmetry in the universe. In addition, near- 
future experiments, such as NO^A, will aim at the de- 
termination of the mass hierarchy [l2j . 

The ongoing experiments might also be able to ver- 
ify observations which do not seem to fit into the three- 
neutrino framework [l3l Il4| . If these revelations were 
confirmed, their explanation would require invoking ster- 
ile neutrinos [15T - [l7l ] and extending the standard model. 

The correct interpretation of the outcome of oscil- 
lation experiments requires a precise knowledge of the 
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(anti)neutrino cross sections. This is the case even in 
those experiments in which the event yields in near and 
far detectors are compared, because, in general, flux and 
backgrounds do not scale in a simple manner [jq . If 
the cross section used in the analysis is underestimated 
in some energy region, a false oscillation signal will be 
obtained as a consequence. 

The description of nuclear effects in neutrino scattering 
is now generally regarded as one of the main sources of 
systematic uncertainties in oscillation experiments. In 
particular, profound understanding of charged-current 
quasielastic (CCQE) interaction with nucleons bound in 
the nucleus plays a crucial role. In the energy region of 
~1 GeV, it is the dominant reaction mechanism. For the 
energy of a few GeV, A-resonance excitation becomes 
equally important llSjl . However, its description involves 
all — but not only [20| — the difficulties appearing in the 
CCQE case. 

The carbon nucleus is of special importance, because 
the targets applied in neutrino detectors often involve 
carbon compounds, such as mineral oil, polystyrene, or 
organic scintillators. In addition, the proper understand- 
ing of carbon may provide useful clues for the description 
of oxygen, as nuclear effects in these two targets are ex- 
pected to be similar, based on the experience gained in 
electron scattering. 

The total 1 gC(^ A1 ,^~) cross section as a function of 
energy has recently been measured with the NOMAD [2l| 
and MiniBooNE (T8J experiments. 

In NOMAD, studying neutrinos of energy down to 
3 GeV, CCQE events with and without knocked-out pro- 
ton detected have been analyzed separately and, after ad- 
justing the description of final-state interactions, shown 
to yield consistent results. Thanks to the 45-GeV average 
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energy of the deep-inelastic-scattering events, the nor- 
malization has been determined from the well-known to- 
tal inclusive charged-current (CC) cross section and from 
the purely leptonic process of inverse muon decay. The 
observed reduction of the cross section due to nuclear 
effects is ~4%. 

In MiniBooNE, all events without pions detected have 
been classified as CCQE and the total cross section is 
extracted for energy up to 2 GeV. The high statistics 
of CCQE events (-10 x those in NOMAD) allowed for 
obtaining, for the first time, the double differential cross 
section. Surprisingly, the reported total cross section for 
carbon is higher than that for free neutrons. The size of 
the effect is —5% at a few hundred MeV, increasing to 
—15% at neutrino energy higher than —850 MeV. 

As a 20% uncertainty of the cross section would have 
an important impact on determination of oscillation pa- 
rameters [IH, the difference between the NOMAD and 
MiniBooNE results requires very careful theoretical anal- 
ysis, in addition to the planned investigation with the 
SciNO^A experiment [231 ]. 

The MiniBooNE Collaboration has also measured the 
cross section for neutral-current elastic (NCE) interac- 
tion, in which nuclear effects are expected not to dif- 
fer from those in the CCQE channel. The NCE data 
may thus provide a cross-check for the applied descrip- 
tion of the target nucleus and, what is even more im- 
portant, additional information on nucleon kinematics in 
the MiniBooNE experiment, unavailable in CCQE scat- 
tering. Therefore, in this paper, we pay special attention 
to NCE scattering. 

A complete description of a nucleus should account 
for both the shell structure and nucleon-nucleon corre- 
lations of various kind. The shell structure is particu- 
larly important for studies involving low-energy probes. 
Correlations, causing strength fragmentation and partial 
depletion of the shells, modify the nuclear structure and 
give rise to deeply bound nucleons with high momentum, 
which are only accessible to high-energy probes. 

Owing to the complicated nature of nuclear interac- 
tions and to the large number of nucleons, in the de- 
scription of the carbon nucleus, approximations are nec- 
essary. One needs, however, to be aware of the accuracy 
of the approximations applied. For example, the nuclear 
transparency of carbon calculated without accounting for 
correlations is underestimated by —15% [241 [25|. 

The relativistic Fermi gas (RFG) model neglects both 
the shell structure and correlations, treating the nucleus 
as a fragment of noninteracting nuclear matter of uniform 
density in a constant potential. Due to its simplicity, it 
is applied in Monte Carlo generators used in data analy- 
sis. One needs to bear in mind, however, that obtaining 
good agreement between the simulation and the distribu- 
tions of collected events may require introducing various 
ad hoc modifications, as a remedy for the shortcomings 
of the RFG approach [H [H, [23. Obviously, all mod- 
ifications applied to experimental analysis must be jus- 
tified on physics grounds within a consistent theoretical 



scheme. 

The MiniBooNE data have been analyzed within var- 
ious models. Martini et al. [28] have developed an ap- 
proach based on the random-phase approximation, which 
subsequently has been improved by introducing relativis- 
tic corrections [2!| . Mainly due to sizable contributions of 
the two-particle-two-hole and three-particle-three-hole 
processes, the results of Ref. [29( are in excellent agree- 
ment with the MiniBooNE data for NCE and CCQE scat- 
tering, including those for the double differential CCQE 
cross section. Analysis of data from other experiments 
in the approach of Ref. [29( is not yet available. It is 
somewhat surprising that for CCQE antineutrino inter- 
actions, Martini et al. predict significantly lower con- 
tributions of multinucleon processes and in the energy 
region of —1 GeV, the total cross section being similar to 
the RFG one. 

This effect is not confirmed by the results of Nieves 
et al. [30I l3~lj , obtained using an effective interaction de- 
termined from data for photon, electron, and pion scat- 
tering off nuclei. Calculations of the double differential 
CCQE cross section have the authors of Ref. [3l| to 
the conclusion that the multinucleon contribution may 
effectively be accounted for by increasing the value of 
the axial mass. It is of great importance, e.g., for Monte 
Carlo simulations, allowing for incorporation of compli- 
cated processes in a simple manner. In addition, com- 
paring their results to the cross section extracted from 
MiniBooNE, Nieves et al. have deduced that the neu- 
trino flux in MiniBooNE seems to be underestimated by 
—9%. As in the calculations of Ref. [3l[ the multinu- 
cleon strength is constrained by precise data, this result 
points to a possible source of discrepancy between the 
cross sections from NOMAD and MiniBooNE. 

Amaro et al. [32j have observed that the relativistic 
mean-field approach provides a very good description 
of the shape of the double differential cross section ex- 
tracted by MiniBooNE, but fails to reproduce its normal- 
ization. The authors of Ref. [32[ have hypothesized that 
accounting for meson-exchange currents could reduce the 
discrepancy. 

Using the spectral function of carbon of Ref. [33| , ap- 
plied also in this paper, Benhar et al. [34j have shown 
that the MiniBooNE data for the differential cross sec- 
tion, be it as a function of muon kinetic energy or its pro- 
duction angle, may be described with excellent accuracy 
when the value of the axial mass is adjusted. However, it 
turned out that the total cross section's dependence on 
neutrino energy fails to reproduce the data. This problem 
has been interpreted in Ref. [34J as likely to be related 
to the flux-unfolding procedure, required to extract the 
total cross section. 

In an effort to understand the MiniBooNE data, 
Meucci et al. [35|, |36[ have shown that the relativistic 
Green's function model may provide a very good descrip- 
tion of the total CCQE cross section while the (double) 
differential NCE (CCQE) cross section is in fairly good 
agreement with the experimental points. The authors 



3 



of Refs. [H, [36| have ascribed the enhancement of the 
cross sections to rescattering of the knocked-out nucleon, 
which may also produce A resonance, and multinucleon 
processes. Within the approach based on the relativistic 
distorted-wave impulse approximation, Meucci et al. [36| 
have obtained the results in somewhat worse agreement 
with the NCE data, but similar to those of Ref. [13] • 

Recently, Lalakulich et al. |38[ have investigated the 
role of multinucleon effects in CC neutrino interactions 
using the Giessen Boltzmann-Uehling-Uhlenbeck event 
generator. In their calculations, the two-nucleon strength 
is fit to the MiniBooNE data using well-motivated phys- 
ically ansatz. Note that in Ref. (38[ the energy spectra 
of the knockout nucleons are presented in detail. 

Th is p aper provides an extension of the approach of 
Ref. [39|, where NCE scattering off oxygen has been 
considered in the context of neutrino-induced gamma- 
ray production, to the carbon target. Some of the re- 
sults for CCQE interaction have been previo usly pre- 
sented (with arbitrary normalization) in Ref. [4U |. dis- 
cussing various issues related to the Q 2 distribution of 
event yield. Here we cover both NCE and CCQE pro- 
cesses, treating them in an identical manner and compar- 
ing the obtained cross sections to those measured by the 
Brookhaven National Laboratory Experiment 734 (BNL 
E734) jU, MiniBooNE [H S3, and NOMAD [H|. 

Our numerical results differ from those presented in 
Ref. [Hj], owing to an error found in the code used to 
obtain the cross sections presented there. 

We acknowledge that the energy dependence of the 
cross sections for neutrino scattering off carbon has been 
also reported by the KARMEN [44] and LSND [H Col- 
laborations. However, in the energy region of these ex- 
periments, ~30-50 MeV and -120-220 MeV for v e and 
data, respectively, scattering off the nucleus is known 
to be dominated by transitions between discrete states 
below the nucleon emission threshold. As we focus on 
nucleon knockout to the continuum, the description of 
such a process lies beyond the scope of this work. 

In Sec. [IT] we consider in the most general manner the 
kinematics of nucleon knockout from a nuclear target, 
recalling the definitions of reconstructed Q 2 applied in 
experimental analysis of NCE and CCQE neutrino inter- 
actions. Our approach is detailed in Sec. IIIII In Sees. HVl 
and [V] our results for NCE and CCQE scattering, re- 
spectively, are compared to the available experimental 
data. In Sec. IVII we discuss how nuclear effects influ- 
ence the determination of the axial mass, and consider 
possible sources of discrepancy between different experi- 
ments. Finally, Sec. I VIll briefly summarizes our findings 
and states the conclusions. 



II. RECONSTRUCTION OF Q 2 IN NEUTRINO 
INTERACTIONS 

In a nucleon knockout from a nucleus at rest by a neu- 
trino of energy E k , the energy and momentum conserva- 



tion read, in general, 

E k + M A = E k , +E A _ 1 + E p , , (1) 
k = k'-p + p', (2) 

respectively, where E p > = \/ M 2 + p' 2 is the energy of 
the knocked-out nucleon of momentum p' . In the case of 
NCE (CCQE) interaction E k < and k' denote the energy 
and momentum of the final neutrino (charged lepton of 
mass m'). The energy of the residual nucleus Ea-i may 
be cast in the form 

E A -i = y/{M A -M + E) 2 +p 2 , 

with M (M A ) being the nucleon (target nucleus) mass. 
The removal energy E is defined as the excitation energy 
of the residual nucleus in its rest frame. The momentum 
of the residual nucleus is denoted as — p, to make use of 
a simple interpretation of p, being the struck nucleon's 
momentum when a neutrino scatters off a single nucleon 
in the nucleus. 

In theoretical neutrino physics, the differential cross 
sections with respect to four-momentum squared Q 2 , 

Q 2 = (p' - p) 2 - [Ep, + E A _! - M A f (3) 
= (k-k>) 2 -(E k -E k ,) 2 . (4) 

are often considered. However, invariant Q 2 cannot be 
calculated from measurable quantities and therefore in 
data analysis reconstructed Q 2 , defined as 

Q 2 CCiNC = 2A/(^-M) (5) 

and 

Q 2 CCiCC = 2E™(E k , - |k'| costf) - m' 2 (6) 

in NCE and CCQE interactions, respectively, is used in- 
stead [3, Hl|, |42j • The reconstructed energy of the in- 
coming neutrino E^, cc , appearing in Eq. ©, is given by 

E roc = 2E k ,M-(m' 2 + M 2 -M 2 ) 
2(M-E k > + |k'|cos0) ' 

where M — M — e, with e being the average value of the 
removal energy. 

The value of the Q 2 reconstructed in NCE (CCQE) 
interactions corresponds to the invariant Q 2 for the 
scattering off a free (bound) nucleon at rest, yielding 
the measured value of the knocked-out nucleon energy 
(charged- lepton energy and its production angle 9) , com- 
pare Eqs. © and © with Eqs. (J)-®. 

Differences between invariant Q 2 and Q 2 ec qq are dis- 
cussed in Ref. .40]. Note that as in Eq. the expres- 
sion [E k r — |k'| cos#) may be greater than M, the recon- 
structed energy does not have to be positive, leading to 
the negative values of Q 2 . cc cc for certain kinematics. 
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III. DESCRIPTION OF THE APPROACH 

We assume that the process of neutrino-nucleus inter- 
action involves a single nucleon and the remaining (A—l) 
nucleons act as a spectator system. This scheme, called 
the impulse approximation (IA), is valid when the mo- 
mentum transfer |q| = |k — k'| is high enough for the 
probe to distinguish individual nucleons in the nucleus, 
as the probe's spacial resolution is ~ l/|q| (see Ref. [Hj]). 

In the IA regime, the neutrino-nucleus cross section 
is obtained convoluting the elementary neutrino-nucleon 
cross section with the hole and particle spectral functions 
(SFs) of the nucleus. 

The hole spectral function Phoio(p, E) is the probabil- 
ity distribution of removing a nucleon of momentum p 
and leaving the residual nucleus with excitation E in its 
rest frame. The particle spectral function P pa rt(p', T 1 ) 
describes the propagation of a nucleon of momentum p' 
and kinetic energy T' . 

In the relativistic Fermi gas model, frequently used in 
Monte Carlo simulations, the hole and particle SFs read 



P£ l(p,E) = 



47174 



(p F - |p|) 5{E p - e - M + E), 



P™ t (p':T')= l-e( PF -\p'\) 8{E P ,-M-T'), 



where E, p — y/ M 2 + p 2 , the average removal energy is 
denoted by s, and pp is the Fermi momentum. 

Note that, in general, the proton and neutron SFs may 
differ. However, for symmetric nuclei (N = Z) it turns 
out that the proton and neutron momentum distributions 
are the same and the difference between the SFs amounts 
to a shift of the removal-energy distributions, 
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taking into account that neutrons are more deeply bound. 
In the case of the carbon nucleus A = 2.76 MeV [47(. 

Realistic hole spectral functions for various nuclei have 
been obtained by the authors of Ref. [33| in the local 
density approximation (LDA), relying on a premise that 
short-range correlations between nucleons are unaffected 
by surface and shell effects. This assumption is sup- 
ported by the observation that for nuclei with mass num- 
ber A > 4, the momentum distribution (per nucleon) 
is largely independent of A for momentum higher than 
~300 MeV [H. In the LDA scheme, the shell struc- 
ture of the nucleus, deduced from experimental (e, e') 
data [49|, |50| . can consistently be combined with the 
correlation contribution obtained from theoretical cal- 
culations for uniform nuclear matter at different densi- 
ties [33|, [5l|. The carbon SF of Ref. [33| has been ex- 
tensively used in the analysis of electron scattering data 
in various kinematical regimes. Moreover, it provides 
a quantitative account of the nucleon momentum dis- 
tributions extracted from (e, e'p) data at large missing 
energy and momentum [251 ] . 



In the IA approach, the total NCE cross section reads 

da NG 

J d»pdEF& le {p,E] 

N—p. n 1 



dQ 2 



£ Jd 3 pdEP» le (P,E)^, 
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where cr^Jy is the elementary NCE neutrino-nucleon cross 
section, 
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with the Fermi constant G F = 1.16637 • ICT 5 /GeV 2 . Al- 
though the Pauli blocking factor is not written explicitly, 
we account for the particle SF using the LDA result of 
Ref. (U- 

The contraction of the leptonic and hadronic tensors, 
L^ v and W^ u , may be cast in the form 

5 

L^W" =2Y,AiWi, (11) 

j'=i 

where Aj are kinematic coefficients equal to 

Ai = 2M 2 k ■ k', 

A 2 = 2p-kp-k' -M 2 k-k', 

A 3 = 2(p-k' k-q-p-kk' -q), (12) 

A A = 2k-qk'-q-k-k'q 2 1 

A5 = 2(p ■ k k! ■ q + p ■ k k ■ q — k ■ k' p ■ q) , 



with p — (E p , p) and q — (E p i — E p , p' - p). The 

re f 
factors by 



structure functions Wj are related to the nucleon form 
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with t = — g 2 /(4Af 2 ). In the case of antineutrino scat- 
tering, the sign of W3 is reversed. The form factors ap- 
pearing in NCE scattering (T^ , Fj[ , and Fp) can be 
expressed by those known from the electromagnetic and 
CCQE interactions [Ff, Fa, and Fp) as 



F l N = ±-(F?-Fp)-2sin 2 w F t N , 



1 As ± g A 



F% = -{F° A ±F A ) = - 



2(\-q 2 lM A y 



(14) 
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(see Refs. |52| - |54| ). In the above equations, the upper 
(lower) sign refers to the proton (neutron) form fac- 
tors, the applied value of the weak mixing angle Ow 
is such that sin 2 6w — 0.23116, m w is the pion mass, 
gA = —1.2673, and the strange quark contribution to the 
axial form factor is set to As = —0.08 55]. As far as 
the strange contributions to are concerned, in our 
approach they are neglected on the basis of experimen- 
tal evidence jH, [13] , on the one hand, and lack of reli- 
able theoretical guidance on their Q 2 dependence, on the 
other hand. 

The results presented in this article are obtained using 
state-of-the-art parametrization of the measurable Sachs 
form factors G E and G% from Refs. [H,^, related to 
the electromagnetic form factors Ff by 
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The dipole parametrization of the axial form factor 
is employed. Calculating the cross sections on free nu- 
cleons, we use the world average axial mass value of 
1.03 GeV [60j . in good agreement with the recent analy- 
sis [6l| of deuteron measurements [62-64]. In scattering 
on nucleons bound in the carbon nucleus, the effective 
Ma = 1-23 GeV is applied, which has been determined 
by the MiniBooNE Collaboration from the first shape 
analysis of the Q 2 ec distribution of the largest statistics 
of CCQE events collected to date [26]. It seems to be well 
justified, as for given Ma the approach of Ref. [26| and 
our calculation yield the shape of the MiniBooNE flux- 
averaged Q 2 CC distribution of CCQE events in very good 
agreement over a broad range of Q 2 ec [io|. While this 
purely phenomenological method of accounting for mult- 
inucleon processes in the carbon nucleus cannot be ex- 
pected to be accurate in any kinematical setups, its valid- 
ity is quantitatively supported by the results of Nieves et 
al. [3l| for the double differential cross section in a broad 
kinematical range of MiniBooNE. In addition, it predicts 
nuclear effects in CCQE and D p scattering to be sim- 
ilar, consistently with the conclusion from the NOMAD 
experiment (2l| . 

As neutrinos may undergo CCQE scattering only off 
neutrons, the cross section reads in this case 
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where the elementary CCQE cross section er cc may be 
obtained from Eqs. (fTUj) - (fll?)) . replacing 



Gf — > Gf cos9c, 



E k , = |k'| -> E k , = Vrn' 2 + k' 2 , 
jf -> Ff - Ft, 



(17) 



F% -+ F A . 



The value of the Cabibbo angle 6c corresponds to 
cos 9 C = 0.97418. 



It is worth noting that the differences between the ex- 
pressions for the NCE cross section given in this paper 
and those in Ref. (65|, describing nuclear effects within 
the relativistic Fermi gas model, can be traced back to 
a different treatment of the off-shell dynamics. 

In our approach, off-shell effects are accounted for by 
using the de Forest approximation [66| . which amounts 
to describing the process of scattering on an off-shell nu- 
cleon as if it involved a free nucleon which absorbed only 
part of the energy transferred by the probe. Note that 
in Eq. (fTTj). even the terms vanishing on the mass shell 
are kept. For example, the structure function W4 con- 
tributes to the NCE cross section only due to off-shell 
effects, because when E p > — E p — Ep. — E k i the coefficient 
A4 is proportional to the neutrino mass. Neglecting W4 
increases the total ^\C{v, v') cross section by 1.4% at 
Ek = 450 MeV. For higher and lower energy the effect, 
resulting mainly from the absence of Fp, is larger, reach- 
ing 3.0% at 100 MeV and 2.0% at 5 GeV. 

The strength of final state interactions (FSI) in neu- 
trino scattering can be deduced from comparison between 
electron scattering data and theoretical results, although 
the picture is obscured by flux averaging. 

The effect of FSI is shown to be sizable in the rela- 
tivistic and semirelativistic approaches of Refs. [67T4691 ] . 
which account for it by means of strong relativistic po- 
tentials or their semirelativistic equivalents. Yielding the 
overwhelmin g co ntribution to the high-w tail of the (e, e') 
cross section [67( , FSI play in these mean-field models an 
essential role in bringing the results in good agreement 
with data [32| . 

Within the SF approach, FSI are much weaker. The 
tail of the cross section originates from short-range cor- 
relations between nucleons in the initial state and is only 
slightly enhanced by FSI (46|. As observed in Ref. [34| . 
in most of the kinematical region spanned by the Mini- 
BooNE data set, the main effect of FSI on the double dif- 
ferential (e, e') cross section is a shift of the quasielastic 
peak of the order of 10 MeV. Using the energy-dependent 
and ^-independent fit to the Dirac optical potential of 
carbon [7(J [7l[ , we find the shift averaged over the Mini- 
BooNE kinematics to be 6.8 MeV, consistent with the 
observation of the authors of Ref. [34|. At higher mo- 
mentum transfers, allowed in NOMAD, FSI are known 
to lead to a redistribution of the inclusive strength, re- 
sulting in the quenching of the quasielastic peak and the 
enhancement of its tails 0, |68| . 

(Quasi)clastic processes constrain the high-|q| contri- 
bution to appear solely at high Q 2 , making it negligibly 
small due to the nucleon form factors. For example, while 
at a neutrino energy of 100 GeV momentum transfers 
up to 100.85 GeV are allowed, as much as 89.8 (97.5)% 
of the 1 |C(t , M , fi~) cross section comes from |q| < 1.645 
(2.747) GeV, the range of |q| for E k = 1 (2) GeV. Hence, 
in the context of our work, the kinematical setup of 
NOMAD does not differ significantly from that of Mini- 
BooNE. 

In this paper, we consider the total cross sections and 
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the flux- averaged differential cross sections da/dQ 2 for 
neutrino and antincutrino (quasi)elastic scattering. As 
FSI may cause only a redistribution and a shift of the 
strength, they do not affect the total cross section. The 
differential cross sections da/dQ 2 are expected to be 
modified by FSI at Q 2 < 0.15 GeV 2 Hi □ because only 
at this kinematics the real part of the optical potential 
significantly changes the typical energy of the knocked- 
out nucleon, compare to Eq. J5]). However, in the low-Q 2 
regime, the validity of the impulse approximation, un- 
derlying our calculations, becomes questionable |40l [72|. 
Therefore, in this article, FSI are not taken into account. 
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IV. COMPARISON TO THE NCE DATA 

From the neutrino experiments which obtained results 
on NCE scattering off nuclei, the highest event statis- 
tics to date have been collected by the BNL E734 and 
MiniBooNE Collaborations, allowing for extraction of the 
differential cross sections. Earlier experiments [73l - l79| . 
aiming at verification of the theory of electroweak inter- 
actions, reported the total NCE to CCQE event numbers 
rates for neutrinos and antineutrinos. 

In the BNL E734 experiment [4lJ, studying vp and 
Dp NCE interactions, the target was composed in 79% 
of protons bound in carbon and aluminium and in 21% 
of free protons. The mean value of neutrino (antineu- 
trino) energy was 1.3 GeV (1.2 GeV). Determination of 
the v and v fluxes involved fitting to the CCQE event 
sample [8(3], under the assumption that the CCQE cross 
sections can be described as in the paper of Llewellyn- 
Smith [8 if with corrections for Fermi motion and Pauli 
blocking |82j. The dipole parametrization of the vector 
and axial form factors was applied, with the cutoff masses 
M v = 0.84 GeV and M A = 1.03 GeV, respectively. 

Based on 1686 (1821) candidate vp (Dp) events surviv- 
ing the cuts, the flux-averaged differential cross sections 
da NG /dQ 2 for scattering off the BNL E734 target were 
extracted [4l| . The normalization uncertainty of the v 
(D) result is estimated to be 11.2% (10.4%). 

We calculate the cross sections for the BNL E734 tar- 
get, neglecting the contribution of aluminium nuclei, ac- 
cording to 



da NC 
dQ 2 



_ j.da up , h _|_ q ^ a v V ,c 



dQ 2 



dQ 2 



(18) 



where / = 0.21 is the fraction of free nucleons in the 
target, whereas a up h (°Vp, c) denotes the cross section 
for scattering off free protons (protons bound in carbon) 
averaged over the neutrino energy 0.2 < Ek < 5 GeV. 

Figure [T] shows that the SF calculations provide a fairly 
good description of the BNL E734 data [41| in the whole 



1 Note that this observation is largely independent of neutrino en- 
ergy and therefore it applies also to flux-averaged results. 



FIG. 1. (Color online) Differential cross section a NC /dQ 2 for 
NCE v and v scattering in the BNL E734 experiment. Our 
results are compared to the data from Ref. [4l]]. The error 
bars do not include the normalization uncertainty of 11.2% 
(10.4%) in the v ip) case. 



range of Q 2 . Note that the agreement seems to be bet- 
ter with the higher-statistics antineutrino data. This is 
also the case for neutrinos in the region of the lowest 
uncertainty, 0.5 < Q 2 < 0.8 GeV 2 . 

We checked that using Q 2 CC defined in Eq. ([5]) instead 
of Q 2 does not change the conclusions significantly, as 
the main difference between these two variables appears 
for the values lower than 0.13 GeV 2 . 

The fluxes, and thus the absolute normalization of the 
BNL E734 data, were determined using a description of 
CCQE interaction clearly different from ours. Do, there- 
fore, our calculations agree with the data in just an ac- 
cidental way? The flux- averaged 1 |C(^, v') cross section 
obtained with Ma = 1-23 GeV varies by a factor of 9.1 
at 0.15 < Q 2 < 1.15 GeV 2 . However, over this interval, 
it differs by less than 10% from the free cross section cal- 
culated with Ma = 1-03 GeV. The corrections for Pauli 
blocking applied by the authors of Ref. [4l[ could dimin- 
ish this difference, which suggests that numerically our 
results may not be very different from those used to find 
the normalization of the BNL E734 cross sections. 

We would like to acknowledge that the BNL E734 data 
have been analyzed by the authors of Refs. [83], H3| to ob- 
tain constraints on the strange contribution to the form 
factors. 

The MiniBooNE experiment, using the Cherenkov de- 
tector filled with mineral oil (CH2), is sensitive to both vp 
and vn NCE scattering [42[. In neutrino mode, its beam 
has been composed almost exclusively of muon neutri- 
nos, with an average energy of 788 MeV. The neutrino 
flux at the detector has been determined by a Monte 
Carlo simulation covered in detail in Ref. [85| • The accu- 
racy of the shape determination has been proven by the 
comparison of the observed and predicted v^ energy dis- 
tribution in the CCQE event sample. The calculations 
were performed within the RFC model with enhanced 
Pauli blocking effect. The dipole parametrization of the 
axial form factor was used with Ma = 1-23 GeV, as in 
this paper. However, the normalization of the measured 
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FIG. 2. (Color online) Differential cross section da NC /dQ 
[Panel (a)] and da^/dQ 2 [Panel (b)] for NCE neutrino scat- 
tering off CEb averaged over the MiniBooNE flux. The spec- 
tral function calculations (dashed line) are compared to the 
data [421 ] with (circles) and without (squares) the correc- 
tion (|22|l . The error bars do not account for the normaliza- 
tion uncertainty of 18.1%. The absolute normalizations agree 
when our results are multiplied by 1.2 (solid line). 



distribution of CCQE events is reported to be higher by 
a factor of 1.21 ± 0.24 than the calculated one 



The MiniBooNE Collaboration recorded in neutrino 
mode 94531 candidate events surviving the NCE cuts, 
which allowed for extraction of the differential cross sec- 
tion with unprecedented precision. 

The theoretical estimate of the differential NCE cross 
section per nucleon as a function of Q 2 reads for the 
MiniBooNE target 



da NC 2 da 
dQ 2 



14 dQ 2 



6 

14 



da 



uN.C 



N—p 1 n 



dQ 2 



(19) 



The most important difference with respect to Eq. (fT5|) 
is the contribution of the neutrons bound in the carbon 
nucleus. 

The detection efficiency for each contributing process 
is estimated in Monte Carlo simulations to be differ- 
ent; therefore the MiniBooNE Collaboration reported in 
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FIG. 3. (Color online) Same as Fig. [5] but on a logarithmic 
scale. 



Rcf. 42] both the cross section, 



da NC 
dQlc 



2 „ da uPi H 

Ti vp ' H ~dW 



6 

14 



da 



uN,C- 



vN, C 



dQ 2 c 



(20) 

and the dependence of the efficiency corrections Cj on 
Q 2 ec , being defined as in Eq. ([5]). Note that in scattering 
off a free proton Q 2 . cc = Q 2 . 

Within a given approach, it is possible to obtain a sin- 
gle effective correction to the theoretical calculations C, 
such that 



dQ 



2 

rcc 



dQ 2 J 



(21) 



with da NC /dQ 2 ec being defined by analogy to Eq. (|T9)) . 

Because the values of the efficiency corrections are 
given only for the values of Q 2 CC corresponding to the 
reported cross section, it is convenient to apply the cor- 
rection to the data rather than to the theoretical result, 
according to 



2 ' 

rcc 



CdQ 



(22) 



In general, such a correction is a model-dependent 
function of the axial mass. It turns out, however, that 
in our calculations, the effects of the proton and neutron 
corrections nearly cancel each other and the value of C 
is very close to 1. The deviations of C from the unit 
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FIG. 4. (Color online) Dependence of the shape of the dif- 
ferential NCE cross section da NC /dQ^ ac in the MiniBooNE 
experiment on the axial mass value. The curves are labeled 
by Ma and the multipliers applied to match their values at 
high Q? cc . 



value, being typically less than 0.25%, exceed 1% in only 
one bin of 0.169 < Q r 2 cc < 0.236 GeV 2 . Therefore, the 
MiniBooNE data with and without the correction (|2"2")l 
largely overlap. 

Figures [5] and [3] show that our calculations reproduce 
the shape of the MiniBooNE NCE v cross section very 
accurately, both in the region of lower Q 2 ec and in the 
tail. For Q 2 ec < 0.64 GeV 2 , the discrepancies are on av- 
erage 1.6%, remaining smaller than 2.7%. The most siz- 
able deviations appear for 0.8 < Q 2 CC < 1.1 GeV 2 , where 
also the data uncertainties are the largest. However, to 
match the absolute scale of the MiniBooNE data, it was 
necessary to multiply the SF results by a factor of 1.2. 
Although this factor is somewhat larger than the data 
normalization uncertainty of 18.1% reported in Ref. [42j . 
it remains consistent with the normalization discrepancy 
of 1.21 ±0.24 observed by the MiniBooNE Collaboration 
in the CCQE analysis [H]. 

Using in comparisons with the MiniBooNE NCE data 
the differential cross section as a function of the theoreti- 
cal Q 2 instead of Q 2 CC , as in Refs. [29|, |H, |86[ , is a fair ap- 
proximation. The differences between them are easiest to 
notice for Q 2 CC < 0.13 GeV 2 , compare panels (a) and (b) 
of Fig. [2J i.e., in the region where no experimental data 
are available. However, when the cross section's slope on 
a logarithmic scale is of interest, the discrepancy, exceed- 
ing 10%, appears also for Q 2 0C > 1 GeV 2 (see Fig.©. We 
observe that for 0.12 < Q 2 CC < 1.75 GeV 2 , the effect of 
using Q 2 QC may be approximated by 



(1 + aQ 2 )- 



dQ 2 



(23) 



with 2% accuracy. When the axial mass of 1.23 GeV 
is applied, the parameter a is equal to 0.095/GeV 2 
(0.070/GeV 2 ) for the SF approach (RFG model). 

It is worth emphasizing that within the approach of 
this paper, it does not seem possible to achieve a good 
agreement with the shape of the MiniBooNE NCE data 
applying an axial mass very different from 1.23 GeV. To 



illustrate this, in Fig. [4] we show the results obtained with 
M A equal to 1.03, 1.23, and 1.39 GeV. Only the middle 
value of the axial mass yields the correct slope of the 
NCE cross section. 

To clarify a possible source of the normalization dis- 
crepancy, it is useful to check how our approach describes 
the CCQE data. 



V. COMPARISON TO THE CCQE DATA 

For CCQE neutrino scattering, the MiniBooNE Col- 
laboration has obtained the differential cross section as 
a function of Q 2 . cc , defined as in Eq. ©, from the 
collected in neutrino mode 146070 events passing the 
cuts [lH . Note that the analysis in Ref. [l8| differs in the 
back gro und treatment from the earlier one, described in 
Ref. [26| (193709 events passing the cuts). 

In Ref. [18[ is also reported the double differential cross 
section da cc jdT^d cos 8, with T M and cos 6* being the 
muon kinetic energy and production angle, respectively. 
As this paper aims at treating NCE and CCQE interac- 
tion on equal footing, we discuss only the data for the 
single differential cross section. 

We have calculated the differential CCQE cross sec- 
tions da cc /dQ 2 and dcr GC /dQ 2 0C , as in Eq. ([15]). and 
averaged them over the MiniBooNE flux. 

Figure [5(a) shows that, as in the NCE case, our result 
needs to be multiplied by 1.2 to match the normalization 
of the MiniBooNE data. The agreement with the Mini- 
BooNE data is fair, but the obtained cross section seems 
to be shifted by —0.05 GeV 2 with respect to the data, 

The shift size coincides 



compare Figs. 5(a) and 6(a) 



with the smallest bin size of the data. 

Our calculations correctly describe the ratio of the 
NCE to CCQE cross section in the MiniBooNE exper- 
iment [42| [see Fig. |5(b)| . Most of the discrepancies are 
related to the shift observed in the CCQE cross section 
[see Fig. |6(b)| . Note that the calculations using theoret- 
ical Q 2 cannot reproduce the shape of the NCE/CCQE 
ratio. 

For the total CCQE cross section, our approach 
yields the energy dependence in good agreement with the 
MiniBooNE result [18| . However, as presented in Fig. [71 
also in this case the factor of 1.2 is required to match the 
normalization. 

The contribution of the A-production events with pion 
undetected, forming an irreducible background to CCQE 
scattering, has been subtracted in the process of extrac- 
tion of the cross section and the associated uncertainties 
have been taken into account in the error bars of the 
data in Fig. [7J For such events, the reconstructed en- 
ergy is lower than the actual value, typically by ^300 
MeV [87| . One may expect their contribution to be 
most pronounced at the energy ~450 MeV, for which 
the cross section receives the background from the peak 
of the MiniBooNE flux. Therefore, in Fig. [7J the differ- 
ence between the upper curve and the central value of 
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FIG. 5. (Color online) (a) The MiniBooNE flux-averaged 
differential cross section for the 1 |C(^ A i,/i~) scattering as a 
function of Q 2 CC (solid line) and Q 2 (dashed line). Our results 
multiplied by 1.2 are compared to the data [lfij. The error 
bars do not include the normalization uncertainty of 10.7%. 
(b) The ratio of the NCE to CCQE cross sections per nucleon 
for CH2 in the MiniBooNE experiment. The line code is the 
same as that for panel (a) . The data are taken from Ref. [12] . 
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FIG. 6. (Color online) Same as Fig.[5]but for the CCQE cross 
section shifted by 0.05 GeV 2 . 
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FIG. 7. (Color online) The total 12 ;C(f M ,/i _ ) cross section. 
Our calculations without modifications (dashed line) and mul- 
tiplied by 1.2 (solid line) are compared to the MiniBooNE 
results , with the error bars showing the total uncertainty. 



the data for low energy is likely to be ascribable to the 
uncertainties related to the subtraction of the irreducible 
background. 

From the analysis presented in this paper, it con- 
sistently emerges that, while correctly describing the 
shape of the MiniBooNE-reported cross sections (be it 
da NC I dQ^ cc , da cc /dQ^ cc , or a cc as a function of neu- 
trino energy), our calculations fail to reproduce their ab- 
solute scale, remaining inadequate by 20%. 

In this context, it is interesting to compare our results 
to the total 1 gC(f M ,^i _ ) and 1 |C(i? AI ,/i + ) cross sections 
measured with the recent NOMAD experiment [2l| . Due 
to the high energy of its broad-band neutrino (antineu- 
trino) beam with the mean value of 25.9 GeV (17.6 GeV), 
the normalization was precisely determined from a large 
sample of the deep-inelastic-scattering and inverse-muon- 
decay events. Using a drift-chamber detector, a total of 
14021 (10358 single track and 3663 double track) neu- 
trino events and 2237 antineutrino events surviving the 
CCQE cuts were recorded. In Monte Carlo simulations, 
nuclear effects were accounted for by using an approach 
based on the RFG model, but with realistic momentum 
distribution applied instead of the step function. 

In the NOMAD experiment, the neutrino events with 
more than two tracks (the charged lepton plus one pro- 
ton) detected would have been removed from the CCQE 
event sample. On the other hand, in the MiniBooNE 
experiment, all events with no pions detected, including 
two (or more) nucleon knockout events, contribute to the 
CCQE cross section. 

Additional nucleons may originate from interactions 
between the struck nucleon and the spectator system 
in the final state, from the nucleon-nucleon correlations 
in the initial state, and from reaction mechanisms other 
than nucleon knockout, such as those involving meson- 
exchange currents. 

The effect of multiproton states due to FSI on the 
CCQE event sample was estimated in Monte Carlo simu- 
lations of the NOMAD experiment to be very small [2l[ . 
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FIG. 8. (Color online) The total CCQE cross section for muon 
neutrino and antineutrino scattering off carbon. The experi- 
mental data reported by the NOMAD Collaboration are taken 
from Ref. [21J | . The error bars show the total uncertainty. 



to some extent, overestimate the cross section, owing to 
inappropriate description of the low-Q 2 contribution. 

From the total CCQE i/ p cross section, the NOMAD 
Collaboration extracted the axial mass 1.05±0.02(stat)± 
0.06(syst) GeV. We estimate that in the SF approach 
that would correspond to 1.17 GeV, due to stronger 
quenching of the cross section. This value has been ob- 
tained without subtracting the correlated strength and, 
as such, gives the lowest Ma required to fit the NOMAD 
data using our approach. 

Analyzing the shape of the Q 2 oc distribution of double- 
track CCQE ^ events for 0.2 < Q 2 ec < 4 GeV 2 , the 
NOMAD Collaboration found M A = 1.07 ± 0.06(stat) ± 
0.07(syst) GeV. Within the quoted uncertainties, this 
value is consistent with our estimate based on the total 
cross section. 



Because the nucleon-nucleon correlations produce in nu- 
clei overwhelmingly proton- neutron pairs [88|, |89(, they 
may have an influence on the track-based measurement 
of the CCQE cross section of neutrinos, but (in the ab- 
sence of FSI) not on that of antineutrinos. 

Analyzing the nucleons with momentum larger than 
300 MeV0 we find that in the energy range covered by the 
NOMAD experiment, scattering off strongly correlated 
nucleon pairs yields ^6% of the inclusive cross sections. 
This number should be considered as an upper bound 
of a difference between the NOMAD- and MiniBooNE- 
reported cross sections that the approach of this paper is 
able to describe. 

The contribution of two- and multinucleon final states 
to the cross section has been thoroughly analyzed by 
Martini et al. [H, Ell and Nieves et al. 0, Hj]. While 
their results do not fully agree, both groups find that 
they provide significant contributions to the CCQE cross 
section. However, the combined analysis of the NOMAD 
and MiniBooNE data in the approaches of Refs. (28l - [3l| 
is not yet available. 

Figure [8] shows that our calculations of the total inclu- 
sive CCQE Vfj, and Sy, cross sections are in good agree- 
ment with the data obtained by the NOMAD Collabora- 
tion [2l| . We want to emphasize that their normalization 
is not scaled by additional factors. Although the SF re- 
sults are higher by -6% than the NOMAD best fit, this 
difference is less that the ~8% (~11%) systematic uncer- 
tainty of neutrino (antineutrino) data. We observe that 
subtraction of the correlated contribution would bring 
the neutrino result in perfect agreement with the experi- 
mental points. However, this may be a pure coincidence, 
as the analogical difference for antineutrinos cannot be 
explained in the same manner. Additionally, the results 
of Ref. (2|| suggest that the approach of this paper may, 



2 For an event to be classified as a double-track one in NOMAD, 
it was required that the proton momentum exceeds 300 MeV. 



VI. DISCUSSION 

With very few exceptions [I(|[9l|]> the axial form factor 
Fa(Q 2 ) is currently parametrized in the dipole form. Its 
value at Q 2 = is known from the neutron beta-decay 
measurements, whereas the dependence on Q 2 is gov- 
erned by the axial mass Ma- 

Because of the sizable contribution of the axial form 
factor to the (anti)neutrino cross section, the axial mass 
may be determined from the total cross sections or from 
shape fit to the event distribution with respect to some 
kinematic variable, e.g., Q 2 . For free nucleons these two 
methods are equivalent, provided the dipole parametriza- 
tion holds true. For bound nucleons, however, this may 
no longer be the case, due to inaccuracies of the ap- 
plied description of nuclear effects. Note that different 
approaches could yield nearly identical (single) differen- 
tial cross sections and, at the same time, different values 
of the total cross section, or vice versa. 

As an illustrative example, consider CCQE (or NCE) 
neutrino scattering off carbon with the MiniBooNE 
beam. At 0.13 < Q 2 < 2.0 GeV 2 , the shapes of the 
flux-averaged cross sections da cc /dQ 2 (or dcr NC / dQ 2 ) 
obtained within the RFG model (p F = 220 MeV, e = 34 
MeV) and the SF approach differ by less than 2.5%, while 
the discrepancy between their absolute values amounts to 
13%. These features result from the fact that the Q 2 de- 
pendence of the cross section is governed mainly by the 
differential cross section on a free nucleon [see Eq. © or 
(|16|)]. whereas its normalization is determined by the dis- 
tribution of strength in the spectral function. Therefore, 
the axial mass determination from the total cross section 
seems to show stronger nuclear-model dependence than 
its extraction from the shape fit. In addition, the proce- 
dure of flux unfolding, necessary to obtain the total cross 
section, introduces additional uncertainties to the result. 

We made use of this observation, showing for CCQE 
interaction that, within the SF approach, the axial mass 
extracted by the MiniBooNE Collaboration from the first 
shape analysis of the Q 2 ec event distribution [26J is in 
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good agreement with the total 1 gC(^ M , [i ) cross section 
measured with the NOMAD experiment [21j . 

To gauge the influence of nuclear effects on the shape 
of the cross section da cc / dQ 2 cc averaged over the Mini- 
BooNE flux, we have calculated it within the SF ap- 
proach and for the RFG model. In the range 0.25 < 
Qrec — 1-0 GeV 2 , used in Ref. [26| for the determination 
of the axial mass within the RFG model, the SF calcula- 
tion differs by less than 5% from that for the RFG model. 
However, at <3 2 oc = 2.0 GeV 2 this difference increases to 
12%. Therefore, the shape of the differential cross section 
as a function of Q 2 CC is affected by nuclear effects more 
strongly than the shape of the differential cross section 
as a function of Q 2 . We have not observed such behavior 
in NCE scattering, applying the definition of Q 2 CC based 
on knocked-out nuclcon kinetic energy [see Eq. ([5])]. 

In Ref. (2(| , the MiniBooNE Collaboration interpreted 
the higher value of the axial mass, extracted from the 
shape of the Q 2 CC distribution of CCQE(-like) events, 
as an effective method of accounting for the nuclear ef- 
fects neglected in the RFG model. 

Within the impulse approximation, nuclear effects tend 
to reduce the total cross section of muon neutrinos 92] 



and slightly increase the slope of the flux-averaged cross 
section da cc /dQ 2 93 . This behavior is easy to under- 
stand: Apart from resonance excitation followed by me- 
son absorption, the only mechanism contributing to the 
CCQE cross section is scattering off a single nucleon. 
The nuclear targets applied in neutrino experiments are 
stable against the emission of nucleons; thus the nucleon 
separation energy is larger than the kinetic energy. For 
an ^1-GeV neutrino interaction, the reduction of the 
phase space due to the binding is therefore more pro- 
nounced than its increase resulting from Fermi motion. 
Moreover, the separation energy increases Q 2 , reducing 
the energy transfer to the nucleon, which results in fur- 
ther quenching of the cross section by the nucleon form 
factors. For high Q 2 , and therefore for low-energy trans- 
fers, the effect of binding is more significant. Therefore, 
the interpretation of the higher value of the effective axial 
mass proposed by the MiniBooNE Collaboration cannot 
be supported within the I A. 

However, the processes which cannot be accounted for 
within the IA framework have recently been shown to 
contribute to CCQE-like scattering, increasing the total 
cross section and decreasing the slope of the Q 2 CC distri- 
bution of events. Nieves et al. [31] have observed that 
these effects may effectively be described using a higher 
value of the axial mass, confirming the interpretation of 
the MiniBooNE Collaboration. 

Nevertheless, the results presented in Ref. 



31] suggest 



that the neutrino flux in the MiniBooNE experiment is 
underestimated by 9.2 ± 3.5% when M A = 1.077 GeV is 
used in the approach based on the local RFG model. Be- 
cause the total cross sections obtained in the RFG model 
are typically higher by ~10% than those in the SF ap- 
proach, the 20% normalization discrepancy that we have 
observed in comparison to the MiniBooNE data seems 
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FIG. 9. (Color online) Comparison of the total ^C^,/! ) 
cross section calculated within our approach to the data from 
the NOMAD [H[ (triangles) and MiniBooNE [3] (squares) 
experiments. The circles show the MiniBooNE data divided 
by 1.2. 



to be in perfect agreement with the finding of Nieves et 
al [3l|. The idea that the MiniBooNE flux might have 
been underestimated is not new 94]. Should it turn out 
to be correct, this would reconcile the MiniBooNE and 
NOMAD data in a simple way (see Fig- EJ) - 

Note that the MiniBooNE flux estimate is based on 
the extrapolation of the cross section for ^ production 
in p-Be scattering measured on a thin target to a target 
35 times thicker. Such calculations are known to involve 
severe difficulties [95$ . 

In various interaction channels, the ratio of the events 
recorded with the MiniBooNE experiment to those pre- 
dicted in the Monte Carlo simulation exceeds unity. It 
is interesting to note that the size of the ratio is similar: 
1.21 ±0.24 in CCQE scattering [11], 1.23 in CC charged- 
pion production [96|, and 1.58 ± 0.05(stat) ± 0.26(syst) 
in CC neutral-pion production [971 ] . The quoted figures 
refer to the simulations with Ma = 1-23 GeV and the 
parameter k, governing the enhancement of the Pauli 
blocking effect, set to 1.019. We checked that our ap- 
proach and the model applied by MiniBooNE yield the 
total CCQE cross section differing by < 2.7% for the neu- 
trino energy higher than 350 MeV. Therefore, we may 
conclude that our 20% normalization discrepancy is in 
good agreement with the MiniBooNE simulations and 
does not seem to be limited to NCE and CCQE interac- 
tions. 

The latter conclusion is supported by the measurement 
of the inclusive CC cross section with the SciBooNE 
experiment [981 ] . using the same neutrino beam as Mini- 
BooNE. The obtained result is larger than its Monte 
Carlo estimate by a factor of 1.12 (k = 1.0) or 1.29 
(k = 1.022), depending on details of the simulation. 

Our calculations underestimate the absolute values of 
the MiniBooNE-reported cross sections, while being in 
a good agreement with their shape and with the NOMAD 
data. 

In the NOMAD analysis, CCQE v events are defined as 
containing at most one proton detected. When the pro- 
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ton's kinetic energy is measured to exceed 47 MeV (mo- 
mentum > 300 MeV), the event is classified as a double- 
track one. The cross section extracted from the single- 
and double-track events, composing 73.9 and 26.1% of 
the collected sample, respectively, is shown not to differ. 
Therefore, only those two- and multinucleon final states 
(2NFS) which involve additional protons of kinetic en- 
ergy lower than 47 MeV and any neutrons contribute to 
the NOMAD result. Moreover, the contributions of the 
2NFS to the single- and double-track samples seem to be 
equal and do not show energy dependence. 

In the MiniBooNE analysis, CCQE events are required 
exclusively not to involve (detected) pions. Hence, the 
cross sections may be increased by a broader class of 
2NFS, involving two or more protons of kinetic energy 
higher than 47 MeV each. 

As a consequence, the 20% discrepancy between our 
calculations and the MiniBooNE data could, in princi- 
ple, be ascribable to more sizable contributions of 2NFS 
than accounted for in our approach by using the effective 
axial mass of 1.23 GeV. Lacking apparent dependence 
on energy, it is, however, constrained by the NOMAD 
results. Note that 

(i) the discrepancy in the NCE channel stems from the 
final-state nucleons with 50 < T < 650 MeV, where 
T is their kinetic energy in total, and appears to the 
same extent over the whole range of T; 

(ii) because the contributions of 2NFS to the 
MiniBooNE-reported NCE and CCQE cross sec- 
tions seem to be equal, the kinematics of the 
knocked-out nucleons in these two cases should not 
differ significantly; and 

(hi) the NOMAD results constrain the missing strength 
to involve at least two protons with T > 94 MeV. 

The discrepancy between our calculations and the 
NCE cross section measured with MiniBooNE is not lim- 
ited to T > 94 MeV, corresponding to Q 2 ec > 0.177 GeV 2 
and remains constant at the interval 50 < T < 650 MeV, 
which suggests that all the contributing 2NFS channels 
are open below T = 50 MeV. These features do not 
seem to be consistent with the hypothesis that 2NFS 
contribute to the MiniBooNE and NOMAD data in a dif- 
ferent manner, but rather point to the flux uncertainty 
in MiniBooNE being higher than reported. 

VII. SUMMARY 

In this paper, we have applied the spectral function ap- 
proach to describe nuclear effects in (anti)ncutrino scat- 
tering off carbon nucleus, treating in a consistent man- 
ner NCE and CCQE interactions. We have considered 
a broad energy range, from a few hundreds of MeV to 
100 GeV. The dipole parametrization of the axial form 
factor with the cutoff mass 1.23 GeV has been used, as 
determined from the shape of the Q 2 CC distribution of 



CCQE events by the MiniBooNE Collaboration in 
Ref. [26j . This effective method of accounting for two- 
and multinucleon final-state contributions to the cross 
section seems to be justified in view of recent results of 
Nieves et al. (3~0 |. 

It has been shown that our approach provides a fairly 
good description of the NCE neutrino and antineutrino 
differential cross sections da NC /dQ 2 measured with the 
BNL E734 experiment. A good agreement has been 
found with the total CCQE and cross sections re- 
ported by the NOMAD Collaboration. 

Our calculations provide very accurate description of 
the shape of the NCE neutrino differential cross section 
da /dQ 2 cc obtained from the MiniBooNE experiment. 
This result confirms that nuclear effects in CCQE and 
NCE neutrino interactions are very similar, because the 
axial mass applied was extracted from the CCQE event 
sample. It turns out, however, that the absolute value of 
the calculated cross section da /dQ 2 cc underestimates 
the MiniBooNE data by 20%. The same discrepancy 
is observed for CCQE interaction in comparisons to 
the differential cross section da cc /dQ 2 cc and the flux- 
unfolded total cross section reported by the MiniBooNE 
Collaboration. 

The difference between the total 1 gC(^ A1 , /i~) cross sec- 
tion from NOMAD (3 < E k < 100 GeV) and MiniBooNE 
(0.4 < Ek < 2.0 GeV) is sometimes attributed to a siz- 
able contribution of two- and multinucleon final states. 
This reasoning is based on the fact that in the latter 
experiment, nucleons knocked out from the nucleus in 
CCQE neutrino scattering have not been detected. 

We have argued, however, that the NCE data provide 
an indication of the kinematics of nucleons in CCQE scat- 
tering, because the MiniBooNE results strongly constrain 
the allowed differences between nuclear effects in NCE 
and CCQE interactions. In the kinematic region used 
to extract the NCE cross section, we find no evidence 
for the contribution of two- and multinucleon final states 
which would not have contributed also to the CCQE cross 
sections reported by the NOMAD Collaboration. There- 
fore, the discrepancy between the results from the Mini- 
BooNE and NOMAD experiments seems more likely to 
be ascribable to underestimated flux uncertainty in the 
MiniBooNE data analysis. 

In this paper we have observed that, although the 
CCQE differential cross section as a function of Q 2 shows 
very weak nuclear-model dependence, this is not the case 
for da / dQ 2 . cc . This feature may be traced back to the 
definition of Q 2 CC applied to CCQE interaction, as does 
not appear in the NCE case. 
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